GLOBAL STABILITY AND LOCAL BIFURCATIONS IN A 
TWO-FLUID MODEL FOR TOKAMAK PLASMA 

D. ZHELYAZOV *, D. HAN-KWAN t, AND J.D.M. RADEMACHER * 

Abstract. Wc study a two-fluid description of high and low temperature components of the 
electron velocity distribution of an idealized tokamak plasma. We reflne previous results on the 
laminar steady-state solution. On the one hand, we prove global stability outside a parameter set of 
possible linear instability. On the other hand, for a large set of parameters, we prove the primary 
instabilities for varying temperature difference stem from the lowest spatial harmonics. We moreover 
show that any codimension-one bifurcation is a supercritical Andronov-Hopf bifurcation, which yields 
stable periodic solutions in the form of traveling waves. In the degenerate case, where the instability 
region in the temperature difference is a point, we prove that the bifurcating periodic orbits form an 
arc of stable periodic solutions. We provide numerical simulations to illustrate and corroborate our 
analysis. These also suggest that the stable periodic orbit, which bifurcated from the steady-state, 
undergoes additional bifurcations. 



1. Introduction. In this paper we analyze the stabihty and prhnary local bi- 
furcations of a laminar steady state in the following two-fluid model for high and 
low temperature in a tokamak fusion plasma near the scrape-off layer. The model 
equations read 



dtp- ^T-d^,p- -E^ -Vp- 
E = -VV, 



1, 



vV'^p- 



fl.l) 



where E 



{E2, ^Eij , > 0, and are posed in the cylindrical domain 



subject to the Dirichlet boundary conditions 



^(0,X2,t) 



F(ii,.T2,i) -0, 

p±(0) , p±(Li,X2,0 



(1.2) 



where p+(xi) 1 - f^, p^,(xi) 

For = and without the Dirichlet boundary conditions on p*, this system 
has been derived in [3] (L\ = L2 was chosen there), to which we refer for details on 
the model origins. Briefly, p^ model miscible phases of 'hot' and 'cold' plasma with 
constant temperatures > > 0, and V the electric potential, driving via the 
drift velocity E y. B /\B\'^ of all charged particles. The addition of viscous terms on the 
one hand allow to model additional physics by adding diffusion or dissipation; on the 
other hand, it changes the system from hyperbolic to parabolic, whose bifurcations are 
easier to analyze. It turns out that v > Q allows for richer destabilization scenarios. 
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In order to relate our results with the hyperbolic system, we include an analysis of 
the case of small > 0. For the benefit of a significant simplification of the analysis, 
we restrict in this paper to the case of equal viscosity for . 

The introduction of viscosity requires additional boundary conditions. The Dirich- 
Ict boundary conditions on are suitable in this context and helpful for our analysis, 
though in other physical contexts these may not be the right choice. Notably, the 
boundary conditions allow for the laminar steady state 

Pss ^ {pts-.P'^s)-: (1-3) 

for which the electric potential and field vanish, and whose relevance for the system 
was noted in [3] for i/ = 0. If > 0, it is in fact the only steady state that is 
independent of X2. In this paper, we present a detailed analysis of its stability and 
bifurcations for v > Q. For moderate viscosity, the equilibrium is unstable in a 
bounded interval [Ari,AT2] of the parameter AT = T+ - T", see Figured! On 
the other hand, the equilibrium is stable for large enough temperature difference, 
AT > AT2 and also for small enough (including negative) temperature difference 
AT < ATi. The parameter AT is relevant in our analysis since it arises in the 
comoving variable X2 ^ X2 — T~t, while T~ is removed. 

For the hyperbolic case = 0, it turns out that ATi = and, for the spatially 
lowest harmonic eigenfunction, 

where i = L2/L1 is the aspect ratio. At Li ~ L2, that is, £ = 1. this is the instability 
region already found in [3]. It turns out that ly > and Li 7^ L2 allows for much richer 
bifurcation scenarios, and it moreover explains the location of the global stability 
threshold AT* = ALi/n'^ as the limiting linear stability threshold for 1/ = via 

lim AT2 = AT,. 

In fact, this is the upper bound of AT for any linear instability. 

One of the original motivations for this study from [3] with Li = L2 was to find 
subcritical bifurcations at ATi,AT2, which would also explain a difference between 
the local instability threshold AT2 (given by spectral stability) and the global stability 
threshold AT* (essentially depending on a Poincare inequality constant). However, it 
turns out that the bifurcations are always supercritical. 

Coming back to the model origins, the sign of AT can be related to the region 
within the tokamak that is modelled by 'good curvature' (negative AT) and 

'bad curvature' regions, which is consistent with the different stability properties for 
positive and negative AT as noted in [5] . The model captures the Electron Temper- 
ature Gradient instability. The modelling and physical relations to L-H transition 
(see [7]) remain to be understood. "Clearly, the model selection criteria, apart from 
the sound physics behind them, should be based on their capability to reproduce key 
experimental facts such as spontaneous L-H transitions, characteristic intermediate 
regimes (such as dithering) , or hysteresis" |B] . 

In this paper, we pursue a mathematical analysis that may serve as a basis to 
investigate further the relations to physical phenomena. The main results may be 
summarized as follows, see also Figures [T] and WA\ 
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Global stability (Theorem \8.S\) . The steady state pss is globally L^-stable for 
AT < and AT > AT*. Global stability for = in a similar region was al- 
ready proven in P] via an explicit Lyapunov functional given by energy conservation. 
In the case of dissipation, improved bounds give the present result with exponential 
decay. Moreover, the global stability threshold is sharp in the sense that it is realized 
as a limiting linear instability threshold in parameter space. 

Local bifurcations ( Theorems \ 5. 1[ \5.!^) . For a large class of parameter configura- 
tions, including T2/T1 < 2\/2 w 2.8, the following holds. At the stability thresholds 
AT,-, j = 1,2, the critical modes are spatially the lowest harmonics, and the system 
undergoes supercritical Andronov-Hopf bifurcations corresponding to periodic travel- 
ling wave bifurcations with velocity uj = 7r(T+ + T~)/L2- Near the bifurcations, the 
reduced dynamics on a center manifold is the generic normal form. For < v <^ 1 
this always holds at AT2, but not necessarily at ATi. 

The local unfolding of the degenerate case ATi = AT2 proves that the two 
branches of periodic orbits are locally connected, and form an arc of stable peri- 
odic solutions. We numerically corroborate that, further away from this degeneracy, 
secondary instabilities occur along the arc. See Figure [Ol and ^ 

In case L2 3> Ti, the primary instabilities can also be higher spatial harmonics, 
even simultaneously. We thus suspect rich dynamics already at onset, but a detailed 
analysis is beyond the scope of this paper. It is also possible, that as AT increases, 
a sequence of destabilization and restabilization occur through different harmonics. 
Roughly speaking, a heuristic interpretation for the model background would be that 
increasing L2 for fixed Ti, AT introduces richer bifurcations from the steady state. 



Fig. 1.1. (a) Schematic illustration of the main case of a primary 1-instability region in the 
stability analysis of the steady-state pss when including viscosity. The global stability threshold AT« 
is larger than the linear stability threshold, even at u = 0. However, in the limit 1/ the lower 
thresholds coincide, and if in addition i — > 00, then also the upper linear thresholds tend to the global 
ones, (b) Sketch of local bifurcation diagram of the steady state u = with supercritical branches of 
stable limit cycle. Solid line represents stable solutions and dashed lines unstable ones. 

This paper is organized as follows. Section [2] contains numerical computations, 
illustrating the results. In 331 we reformulate the problem for a subsequent bifurcation 
analysis. Section|3]concerns the spectrum of the linearized operator around the steady 
state Pss ■ In fj5] we discuss the center manifold reduction, reduced vector fields and 
prove the main bifurcation results. In ^Sl we explain the relation to travelling wave 
bifurcations, and briefly consider pattern formation in case of an infinite strip. In fj7] 
we discuss nonlinear instability for ly > in the linearly unstable region. Finally, fJS] 
contains the global stability result. 
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2. Numerical results. For illustration of the upcoming analytical bifurcation 
results, we present in this section some numerical computations. We compute the 
deviation u = p ~ pss (see and discrctize with a finite-dimensional spectral 

decomposition (see (|4.3p for the definition of the harmonics gk)'- 

ui{x,t)=Y, C'fci,fc2,i(i)3fci,fc.(x) , /-1,2. (2.1) 

ki = l k2 = -N^^ 

We integrate the resulting system of ODEs using a semi-implicit Crank-Nicolson 
scheme, where only the linear part is implicifl We used the parameter values 

= 9.10"^ ii = La = 2, = 10"\ (2.2) 

while AT = r+ — varied across the instability region. Note that this lies in 
the region L2 < 2V2ii, thus the primary instabilities come from the lowest spatial 
harmonics as proven in Theorem 14.41 All the simulations were made with N^-^ = 
= 32, though we selectively checked with N^^ = N^^ = 64. 




Fig. 2.1. Contour plots of ui{t2) = pi{t2) — Pss with t2 sufficiently large for AT & (ATi , AT2) 
and parameters as in 112.211 . (a) the dynamics is a translation in the periodic X2-direction, AT = 
0.159291, and (b) AT = 0.146122, the dynamics resembles a modulated travelling wave. 

In Figure [^m we plot two periodic travelling wave solutions near the upper stabil- 
ity threshold ATa ~ 0.162 and further inside the nonlinear regime as can be seen by 
the locus of parameters in Figure [221 a)- The weakly nonlinear solution for AT sa ATa 
closely resembles the unstable eigenfunction, while the solution further inside the non- 
linear regime has a clear nonlinear structure. 

In order to trace the stable branches of solutions bifurcating from the supercritical 
Andronov-Hopf bifurcations at AT = ATi, AT2, we perform a simple continuation: 



^Wc modified a code by Jean-Christophc Nave - MIT Department of Mathematics, 
jcnavc@mit.cdu for Navicr-Stokcs equations in vorticity formulation. 
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for AT near the bifurcation at ATi, we simulate an initial condition close to pss 
and after a long transient compute the sup-norm over a long time interval. We then 
slightly increase AT and repeat this step with the initial condition being the solution 
at the final time of the previous step. In this way we obtain the bifurcation diagram 
in Figure 12. 2[ where the numerical instability thresholds are in very good agreement 
with the analytical ones. 




0.05 0.1 0.15 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

AT AT 



(a) (b) 

Fig. 2.2. Bifurcation diagrams in AT with maxtj<t<t2 ||'"i(i)||oo on the vertical axis, where 
ti,t2 are taken sufficiently large (ti > 400J, so that we get a good approximation of the attractor. 
(a) parameters as ii i2.2\ and (b) L\ = 0.9, L2 = 0.9, v = 9.10~3, the diagram is in agreement with 
corollarv \5.S\ 



As predicted by Theorem l5.ll the slope of the resulting curve is larger near the left 
endpoint of the instability region than near the right endpoint. Further away from 
these endpoints, the solution shapes change and we conjecture a period doubling 
bifurcation near AT = 0.0646, rapidly followed by a torus bifurcation. This would 
be consistent with the sharper increase in the sup-norms in Figure 12.21 and the fact 
that the periods of the solutions become rather large, see Figure I^^T b). Despite this 
detour to another attractor, the solution is eventually turning into the near harmonic 
periodic solution bifurcating from the right endpoint AT2, see Figure [^?T] (a). 
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Fig. 2.3. Time traces of ui at the midpoint {xi,X2) = {Li/2, L2/2). (a) AT = 0.056, and (b) 
AT = 0.08. Compare FiaureW^ 
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As an example for instabilities caused by higher spatial harmonics, we plot in 
Figure [231 a solution that emerged from an instability with wavenumber k2 =2. 




0.5 1 



Fig. 2.4. Contour plot of u\(t2) inside the nonlinear regime. In this case, the critical eigen- 
f unction has a wavenumber k2 = 2. 

3. Reformulation and setting. For the bifurcation study it is convenient to 
formulate p.ip through the deviation u = (mi,U2) from p^s, 

p+ = Ul + p+ , =U2+ p^,. 

In terms of u, and in the comoving variable X2 — > X2 + T^t, system reads 
' dtui ^ ATd^^ui + E2/L1 - ■ Vui + i/V^ui, 

< E = -VV, (3.1) 
-VV = U1+U2, 

X e [0,Li] X R/L2I .t > 0, 

subject to (periodic b.c. in X2 and) homogeneous Dirichlet boundary conditions 

Ml(0, X2,t) ^ U2{Q, X2,t)^ 1/(0, X2,t)^ 0, 2^ 
Ui{Li,X2,t) = U2{Li,X2,t) = V{Li,X2,t) = 0. 

Remark 1. We want to briefly point out a peculiarity of the nonlinearity in 
(|l.ip and equivalently (j3.ip .' viewed on complexified phase space, each eigenspace of 
the laplacian is flow invariant and the dynamics is purely linear. 

Indeed, take an eigenfunction e with eigenvalue A and set Uj = ajc with aj € C 
so that E = {oil + Q;2)/AVe. Hence, E^ ■ ^Uj ~ so that (|3.ip is in fact linear. 

However, this does not provide flow invariant spaces for the real equations since 
all eigenvalues and eigenspaces are complex, and the previous argument is incorrect 
for linear combinations. Co-moving frames do not generate real eigenspaces due to 
the asymmetric advection terms. 

Next we choose a simple functional analytic setting for a formulation of (j3.ip as a 
parabolic problem by solving the Poisson equation. This is convenient for the center 
manifold reduction, but also gives a simple well-posedness setting. 
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Let n ■= [0, Li] X [0, L2] and denote the Sobolev spaces = H^{[0, Li] x R/L2Z) 
as well as 



X ■.^H^{[0,Li] X M/L2Z), 

Y := {/ e H2 : /(O, X2) = /(Li, X2) = 0}, (3.3) 
Z:={/eH3 : /(0,a;2) = /(Li,X2)=0}, 

which incorporate the Dirichlet boundary conditions. Wc shall use standard nota- 
tion: for /i, /2 € L^([0,Li] X M./L2Z) wc denote the scalar product by (/i,/2) = 
/n h{x)h{x)dx and for f^ = (fj,i,fj,2) e i'([0,Li] x M/L2^)' j = 1,2 by (fi,f2)2 = 

(fia,f2,l> + (fl,2,f2,2>- 

Thanks to these boundary conditions, we can solve the Poisson equation in (|3.ip : 
see also 21 for explicit solutions. We thus obtain E via the bounded operators Aj : 
Z ^ H3 defined by 

A,/:=a,^(V2)-i/, J = 1,2, 

A/ = (Ai/,A2/f, (3.4) 

A^/=(A2/,-Ai/)^. 

Notably, A2 in fact maps into Z, because E2 — dx2V vanishes for xi — 0, Li due to 
the Dirichlet boundary conditions. 

In order to apply standard results on parabolic equations, let us write (|3.ip equiv- 
alently in the standard form 

fill 

^=Lu + R{u), (3.5) 
at 

so that solutions of this and (jl.ip are in 1-to-l correspondence. Here 

^ fATdx.ui + ^A2iui + U2) + i^V^uA 

"^V -^A2{ui+U2) + l'y^U2) ' 

_ f-A^{ui+U2) ■S/ui\ 
^^''>-[-A^{u,+U2)-Wu2)- 

Note that L G C{Z x Z,X x X) is the linearization of in pss- We have that 
R: ZxZ^YxY since Wuj vanish at xi = 0, Li and is a Banach algebra; 
R is in fact analytic in u. See also ^ Moreover, the imbeddings ^ X^ 
are dense and the uniformly elliptic operator — L : ZxZcXxX^XxX 
is a sectorial operator, generating an analytic semigroup, and so (|3.ip admits mild 
and classical solutions u{t) for any initial condition u{0) £ Y x Y. The sectoriality 
is a consequence of the fact that the laplacian is sectorial in Y with domain of 
the cylinder [5], and this is robust under addition of the lower order terms in L. It 
thus also possesses a square root, which then provides an isomorphism from to 
X. Hence, L is also sectorial on Z with domain X. Note also that L has a compact 
resolvent and thus discrete spectrum accumulating at —00. We discuss its spectrum 
in detail in the next section. 

4. Spectrum of the linearization. For the bifurcation analysis, we distinguish 
the stable spectrum of L, cr-(L) := {A e cr(L) : < 0}, its neutral spectrum cro(L) = 
{A G cr(L) : SRA = 0} and its unstable spectrum cr+(L) = {A e cr(L) : SRA > 0}. 
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The next Lemma characterizes the spectrum and is the basis for the identifica- 
tion of bifurcations. While this concerns the comoving variable of system p.l|) . the 
spectrum for the original system is the same up to a scaling of the imaginary parts. 
See fH 

Lemma 4.1. The spectrum cr(L) o/ L consists of the eigenvalues 



>^k = 17^4^ - ^'t2 ( + ^ ) ± V^, fc e N. X Z, (4.1) 



In 



where N* = N \ {0} and 

particular, G (t_(L), and < then G tT_(L). Moreover, 3? (^A^^ < 

We will start to discuss the relevance and implications of this result after the 
proof. In preparation of the proof, choose the orthogonal basis of X given by 

gk{x) := sm — e ^2 , (4.3) 



where G x Z. In order to express the operator A, denote 

(j)k{x) := COB [— — le ^2 . (4.4) 

Indeed, if / G X, the explicit solution to the Poisson equation — V^T^ = / in terms of 
this basis reads 

We therefore get the explicit formula for A: 

Proof. [Lemma 14. Ij Consider functions of the form qgk{x), fc G N* x Z, where 
g G is an arbitrary constant vector. Since 

A2gk{x) = 3—^-^3,. (x), (4.6) 

the action of L on such functions is 

{Lqgk){x) = Mkqgkix), (4.7) 

where 

/Ci(fc)Ar - C2{k) - C3(fc) -C2(fc) \ , . 

''''=■"1 C2(fc) C2(fc)-C3(fc)j ' ^^-^^ 



with 



2i 



The eigenvalues of are readily computed to be A^. The claims on the real 
parts of immediately follow from inspecting (j4.ip - in particular monotonically 
decreases in fci. □ 

Note that the proof also implies that eigenfunctions of L have the form 



Ck{x) := £,kgk{x) ^ Z X Z, 



(4.9) 



with ^fc a eigenvector of . 

The last statement in Lemma 14.11 means that only Dk > and A^ ^.^^ with 
^2 £ Z \ {0} allow for destabilization, and the real part in this case is given by 



5r(a+ 




(4.10) 



Note that this is a function of the three parameters v/L\, AT/ii, (kij^Y'- ex- 
pected, increasing viscosity always stabilizes, with increasing impact for increasing 
{1^2/ f-Y- However, the dependence of the real part on k2/t is not necessarily mono- 
tone, which allows for intricate destabilization scenarios. 

The imaginary part, 9(A^ ^2))' ^® never zero, which means that all bifurcations 
arc non-stationary and we generically expect Andronov-Hopf bifurcations, where ^2 
determines the wavenumbcr of bifurcating solutions. 

We consider the temperature difference AT as the primary bifurcation parameter 
and therefore focus on the location of instabilities as AT varies, as well as on the 
wavenumber of destabilizing modes determined by k2- 




0.10 



AT 




0.02 0.04 0.06 0.08 0.10 



AT 



(b) 



Fig. 4.1. Real parts of eigenvalues as functions of AT for L\ - 
eigenvalues with k S {1, 10} X {—10, 10} are plotted for (a) y = 10" 



: L2 = 1. The unstable 
3_ (h) V = A- 10^3. 



In Figure WA\ we plot sample computations of spectrum as AT varies, illustrating 
the stabilizing effect of the viscosity. Crossings of eigenvalue curves at zero real part 
can occur, which is expected to generate rich bifurcations. However, in this paper we 
focus on simple Andronov-Hopf bifurcations. 
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Recall the spectral conditions at a primary Andronov-Hopf bifurcation 

(i) There is a constant 7 > s.t.sup{5RA : A € cr_(L)} < —7, 

(ii) cro(L) = {±iuj},uj > and ± iw are simple eigenvalues, (4-11) 

(iii) a+(L)==0, 

and in the nondcgenerate case, the critical eigenvalues transversely cross the imaginary 
axis upon parameter variation. 

It turns out that we can characterize a large part of parameter space, where 
critical eigenvalues have /c2 — 1, that is, k ~ kc := (1,1)- We therefore define the 
following particular case of (|4.11[) . 

Hypothesis 1. It holds that = ~ and there is j > such that 

^X^ < for k en^ xl.\{kc,kc}. 

Here and in the following we denote k = — K2) for k S M^. 

Rearranging sign conditions on (j4.10p and squaring, we readily compute that the 
sign of 9?(A^ for K2 = fci is the sign of 

which is somewhat simpler to handle. In particular, zeros of d are the critical eigen- 
values for bifurcations. This yields the following a priori bounds on AT for linear 
instability. 

Lemma 4.2. For all and K2 > 0, the real roots of d{-,K2) lie in [0,4Li/7r^]. 
Moreover, the real roots approach the endpoints in the limit £ — > 00 if v = o{l~^). 

Proof. Since d(0, K2) < and c)atc?(0, K2) > the lower bound holds. For 
the upper bound, observe that d(4Li/7r^, K2) < and dATd{4:Li/Tr^ , K2) < 0, which 
proves the claim since the quadratic coefficient of AT is negative. The statement on 
the limits readily follows from (|4.12[) upon multiplication by □ 

Note that d{-, K2), as a quadratic polynomial in AT, has two real roots ATi(k2) < 
AT2(k2) if and only if the viscosity is sufficiently small, 

with a double root at equality. Hence, this is a necessary and sufficient condition for 
the occurrence of critical eigenvalues A^i , — \ as AT varies. However, it is subtle to 
determine when the critical eigenvalues destabilize the equilibrium as this requires to 
exclude unstable eigenvalues for all other k2- 

Nevertheless, the location of these parabola's maxima in AT is at 

AT = '^-^ (4 14) 

which is strictly decreasing in K2- Therefore, the fc2-value of these parabola in AT 
can be identified by the relative location of their maxima. 

Remark 2. For H2 = I the roots satisfy ATi = 0{v'^) and AT2 = jfp^^^^^ + 
0{v'^), which was already illustrated in Figure[l[ 
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The geometric nature of bifurcating solutions is determined by the A:2-value of 
critical and destabilizing eigenvalues as AT in- or decreases from outside [0, ALi/tt'^]. 
We thus define 

Definition 4.3. For given Li,£, v, we say that L possesses a k2-instability 
region, ifd{-,k2) has two positive roots Ari(fc|) < AT2(fc|). We call a k2-instability 
region locally primary, if there is a neighbourhood S* C R o/ J := (ATi(fc|), AT2(fc|)), 
s.t. the steady state u ~ is stable for AT G S \ J and ATj(fc|) ^ ATj{K2) for 
K2 7^ fc|, j ~ 1,2. Moreover, we say that the k2-instability region is primary, if it is 
locally primary and S = R. 

To ease notation, we simply write ATj for ATj{l), j ~ 1,2. 

As a first step to understand the nature of destabilizing fc2-instability regions, we 
consider the case ^2 = 1 and in preparation define the following condition. 
Hypothesis 2. Suppose that for given Li,i,v > Q we have 



16^4^3;^, 

and V < i^crrf(l) for AT G {ATi, AT2}, and all ^2 G N*, k2 > 2. 

Note that Hypothesis [2] requires a ratio of temperature difference and viscosity to 
dominate a ratio involving domain geometry and linear mode harmonics. 

Theorem 4.4. 

1. A 1-instability region o/L is locally primary if and only if Hypothesis\^ holds. 
If it holds, then Hypothesis[l\ is satisfied at AT = ATj, j = 1,2. The critical 
eigenvalues are Xj = iiwj with ujj = TTATj/{£Li). 

2. ForO<£< 2V2 w 2.8 any 1-instability region ofL is primary and Hypothe- 
sis{J\ is satisfied at AT = ATj, j = 1,2. 

The point of the theorem is that it provides conditions (Hypothesis [5] or £ < 
2\/2) under which the destabilizing mode for increasing and decreasing AT is known, 
namely the lowest spatial harmonic. Note that the values of £ in particular include 
the case £ — 1 considered in . 

Proof. 

1. A direct calculation gives 

d(AT, 1) - d(AT, K2) 1&ATk2£^LI - (4 + £'^){'iK2 + £'^){£'^ - IQk2)v'^tt^ 



K2-I K2^4(4 + £2)(4^2+^2) 

(4.16) 

In particular. Hypothesis [2] is indeed equivalent to a 1-instability region being 
locally primary. The claims on AT follow readily from inspection of the zeros 
of d. 

2. This is the trivial observation that the right hand side in condition 14.151 is 
strictly negative for these values of £, while the left hand side is positive at 
all possible real roots AT of c?(-, K2) on account of Lemma [ 



□ 

fact 



Remark 3. The critical frequencies in the original X2-variable of (jl.ip are in 



r+ + T^ 
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Now, we are going to present a condition, which guarantees that other destabi- 
lization scenarios also occur. 

Corollary 4.5. Let K2 > 1 and let £ be the unique positive solution £ = of 



K2£^ - 80^2^^ - 64k2(2 + K2) = 0. 



(4.17) 



Then for v = !/crit(l) the 1-instability region is a point, ATi = AT2, that coincides 
with AT2(k2). Notably, £^2 strictly increasing in K2- 

This means that the 1-instabihty region is not primary. In fact, it is also not 
primary for nearby parameter values that produce AT2(4) > AT2(1). The solution 
to for K2 = 4 is £4 « 5.37, and for K2 = 9 it is £g « 7.22. See Figure O 

For £ between these value (and slightly above £g), we numerically find that the 2- 
instability region is primary. We omit the tedious analysis. In general, for any given 
^2 Hypothesis [2] is violated for sufficiently large £ (with j^, Li fixed), since AT is 
bounded (Lemma 14. 2[) . 

Remark 4. It is possible to show that for v small enough, there is a primary 
1-instability region, if £ < £* sa 4.053, where £* is the unique positive root of the 
polynomial 16{£^ + 4)^ - {£^ - 8){£^ + 8){£^ + 16) = 0. 





(a) 



(b) 



Fig. 4.2. Real parts of eigenvalues A^, k = (1,^2), as functions of AT. The parabola are 
ordered in k2 by the decreasing location of maxima. Here v = i'crit(l) that ATi = AT2. (a) 
I = £4, with near primary 2-instability region, (b) I = Ig, with primary 2-instability region. 



corresponding critical value of AT = .^^^jj^p-^ at the double root into the nominator 



Proof. [Corollary 14. 5| Substituting the critical v'^ — ^^^'^j"^i^ti from (|4.13p and the 
ue of AT 

of the right hand side of 14. 161 gives 

r(64K2(2 + K2) + 80k2^^ + 4^2^'' 



4£6Lf 



(4 + ^2)3,^2' 



1, 



where k;2 = ^2- The first factor is positive and roots of the second factor, which we 
denote by g, precisely solve (|4.17p . We have 

d^p^q = 80k2 +8k2^2 



3r 



which is positive at ^ = so that the cubic q with negative cubic coefficient has 
a unique positive root. In addition, this implies that diq < at this root so that 
together with 

d^^^q ^ 4(32 + 32^2 + 20£2 + > 
we infer from implicit differentiation that the location of this root strictly increases 

with K2. □ 
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For the case of small viscosity (and i > 2.8), we omit the somewhat tedious 
detailed analysis for the destabilizing left endpoint. However, we immediately obtain 
the following. 

Corollary 4.6. As — > 0, L has k2 -instability regions for ^2 — > oo with 
ATi(fc|) < AT2(fc|). For sufficiently small v, the conditions of Hypothesis [l\ are 
satisfied at AT2(1), and this is an instability threshold. 

Proof. The presence of all fc2-instability regions clearly holds at = in view of 
(|4.13p . In addition, from (|4.16p we infer at = that 

d(AT,l) - d(AT,K2) > 0, 

so that the critical eigenfunction at the right endpoint of the instability interval has 
mode number ^2 = ±1. This persists for sufficiently small > 0, since the thresholds 
depend continuously on v, and again from 14.161 we see that for each v > there is 
only a finite range of K2 values, for which d{AT, 1) — d{AT, K2) < is possible. □ 

Lastly, we point out the possibility of multiple disjoint primary fc2-instability re- 
gions, where changing AT destabilizes and stabilizes multiple times. In Figure we 
plot eigenvalue curves, where two fc2-iiistability regions consist of a point. Parameters 
I' = i'crit(fc2) = '^01-11(^2) ^-iid i that produce such scenarios can be readily computed 
from (j4.13|) : here we take /c2 = 1, fcj = 4. For perturbed i' < i'crit(l) the instability 
regions become disjoint open intervals. 

Re A 
0.000 

-0.002 

-0.004 

-0.006 

-0.008 

-0.010 

Fig. 4.3. Real parts of eigenvalues A^, k = (1, fc2), as functions of AT. Here £ = 2\/2 + 3\/2 
so that i^crit(l) = i^critC^). In the region plotted only the curves with A;2 = 1,4 are present. 

Remark 5. On account of ()4.13p . for decreasing v and also for increasing £, there 
is an increasingly long sequence of secondary instabilities of Andronov-Hopf type as 
AT increases from zero, with higher and higher spatial harmonics, and another reverse 
sequence as AT reaches AT* . See Figure \41\ 

5. Center manifold reduction. In this section, we consider the vicinity of 
parameters with critical AT = ATj for j = 1 or j = 2 and assume that no other 
eigenvalues lies on the imaginary axis. The main example is a primary 1-instability 
region. For the unfolding of the bifurcation in the generic case ATi < AT2 we 
introduce the parameter /ii by AT = AT^ + fti. In the degenerate case ATi = AT2, 
where ly = t'crit(l)j we additionally unfold with ^2 defined hy ly = i'crit(l) ~ mI- For 
readability we frequently suppress the index j. 

At bifurcation, the critical eigenvalues are then ±ia; and we denote the associated 
eigenfunctions by C{x) := Ck^{x), C(a^): see (|4.9p . Then L possesses a two-dimensional 
real central subspace Ec := spanjSRC? ^^C} and we will show that there is a 

locally invariant 2D center manifold 




Wc = {uq + tp{uQ,fL) : uq G OeJ C 0^2 , G Or2, 
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with V : Oe, El, E^® El = Z^, and neighbourhoods Or2 of /i = 0, and Ob^, 
Oz'i oi G Z^. In case of a primary bifurcation the center manifold is also locally 
exponentially attracting. 

Since we consider k = kc = (1, 1), it is not surprising that the coefhcients Cm{kc) 
defined in (j4.8p show up. It turns out that following modifications are convenient. 



ttAT 2 2/^1-1 , ^ 



We first consider the generic case of (|4.13p . where the unfolding goes by /ii only. 

Theorem 5.1. Assume that Hypothesis]^ holds for a fixed parameter set for 
which ATi < Ar2- Then the steady state u = of system (|3.5p possesses a locally 
exponentially attracting and locally invariant 2D center manifold near u = with the 
reduced dynamics 

dr 

^=iuJz + ^l^az + bz\z\' + 0{\z\{\^i\ + \z\'Y), (5.2) 

where 

ttAT 

_ 27rc^ 
"~ L2(ci-C3i)(cii-2c2i-C3)' (5.3) 
Ll c\ + C3 



h = - 



1^1 



The following corollary proves the nature of the resulting bifurcations, sec Fig. 
[TJa) for an illustration. 

Corollary 5.2. Assume the Hypotheses of Theorem \5. li Then the steady state 
u ~ of system p.Sp undergoes a generic supercritical Andronov-Hopf bifurcations as 
fii varies. Specifically, the reduced vector field coefficients satisfy b < 0, 3(a) = 0{v'^), 
and sgn(5Ra) = -(-1)-' at AT = ATj. 

In particular, near the stability thresholds there exist heteroclinic connections be- 
tween the unstable steady-state and the stable limit cycle. 

As V 0, the radius of the limit cycles, \z{t)\, scales near ATi as \z(t)\ oc 
i^^VAT- ATi, and near AT2 as \z{t)\ (x AT2 - AT. 

Before giving the proof, we formulate the result for unfolding the codimension-2 
case ATi = AT2 , where the critical eigenvalues do not transversely cross the imaginary 
axis. 

Theorem 5.3. Assume that Hypothesis]^ holds for a fixed parameter set for 
which ATi = AT2 ■ Then the steady state u = of system (|3.5p possesses a locally 
exponentially attracting and locally invariant 2D center manifold near u = with the 
reduced dynamics 

dz 

— =\{uj + aoiiijz + ai/^i (02^2 - 03/^1) z + bz\z\^ + TZ (5.4) 
dt 

n^o{^ll + Vl,^Jl\^\z\{\^Ji\^\z\'Y), 
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where aj G M, j = 0, 1, 2, 3, are given by og = 3(a), 

3?(a)7r 1 TT 

L2(C2-Ci) ^/ttLiL2 

and a, b are the constants from Theorem \5. 1\ 

In particular, /or < \fi2 \ ^ I7 there exists a branch of stable periodic orbits, that 
is parametrized by fii and that terminates in supercritical Andronov-Hopf bifurcations 
at AT + ni= ATj, j = 1,2. 

The following Theorem shows, that the bifurcation results 15.1115.31 can be gener- 
alized to instabilities caused by higher spatial harmonics. 

Theorem 5.4. Assume (|4.1ip holds with critical wavenumber k2 , so thatX^^ = 
iu). //Ari(fc|) < Ar2(^2) then the statements of Theorem \5.1\ and Corollarv \5.2\ hold 
with ATj replaced by ATj(fc|), and L2 replaced by iv2/fe for the coefficients on the 
center manifold. If ATi(A:|) = AT2(A:f ) then the statement of Theorem \5.S[ holds for 
the same modifications. 

Proof. Under condition (|4.1ip the center manifold theorem applies as in the first 
parts of the proofs of Theorems 15.11 and 15.31 This yields a stable locally invariant 
manifold with reduced dynamics of Hopf normal form. The only remaining question 
is the sign of the coefficients. 

If k2 is the critical wavenumber in a;2-direction on the domain [0, Li] x [0, L2] then 
1 is this wavenumber on the domain [0,Li] x [0,272/^2] so that Hypothesis [1] holds 
there. Hence, on this domain and with the modifications in the claim. Theorem 15.11 
Corollary 15.21 and Theorem 15.31 hold fully. 

The theorem now follows since the bifurcating branches imbed into the original 
domain. □ 

Remark 6. Recall that there is a sequence of secondary Andronov-Hopf insta- 
bilities as noted in Remark\^ Whenever these occur with a simple pair of complex 
conjugate eigenvalues, analogous center manifold reduction results hold for an un- 
stable 2D manifold. The reduced vector fields are of the same form with coefficients 
given analogous to the above results, but to be computed at different k2 and other 
parameters. 

We start with the proof of Theorem 15.11 

Proof. [Theorem 15.1] For the unfolding with /ii, we modify the definition of R 
in (|3.5p by adding the term ^18x2^1 in the first component and denote the result by 
R{u; pLi). For the resulting bifurcation problem, we verify the hypotheses of the center 
manifold theorem [H Theorem 3.3, p. 46]. 

As noted after (j3.5p . L € £(Z^,X^) is sectorial so that Hypothesis 2.7 in that 
theorem holds, using [H Remark 2.18 p. 37]. Hypotheses 3.1(i) and 2.4 hold on account 
of Theorem 14.41 It remains to show Hypotheses 3.1(ii): smoothness of R. From (|3.5p 
we explicitly compute 

n??^,-„ \„ (liidx2Vi- A^{vi+V2) -^ui- A^{ui+U2) ■yvi\ 

^1)'' = -^AHvi + V2) ■ VU2 - AHu, + U2) ■VV2 )' 

n2o/ ^^ 1 fA-^{vi+V2)-Wwi+A^{wi+W2)-'^Vi\ 

D R{u; ^^,)[v, w] = ~ ^^^^^^^ ^ ^ ^^^^^ ^ . 

Note that A-^{vi + W2) G and Ww E Y^. Since is a Banach algebra (see for 
instance [H Theorem (4.39]), there is a constant Cq > 0, such that 

\\A^{vi+v2)-\/w\\y < Co\\A^ {Vl + i;2)||y2||Vw||y2. 
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Hence \\D^R{u; hi)[v,w]\\y2 < C\\v\\z2\\w\\z2 , that is, Riu; m) G C^{Z^,Y^). More- 
over all the higher derivatives are identically 0, hence R is analytic. This establishes 
the existence of the 2D center manifold and smoothness of as needed below, and for 
which the reduced dynamics has the normal form (j5.2p . Here the critical frequency is 
uj = Ti/ST / L2 due to Theorem 14.41 In order to analyze the coefficients of the reduced 
equation, we write functions in the central subspace as 



Mo(t) = z{t)C + z{t)C, z{t) e C. 

Using the expressions in l4j p. 125] (see also Scholarpedia on Andronov-Hopf bifurca- 
tion), we have 

a = (i?ii(C) + 2i?2o(C, C)2, (5.5) 

b = (2i?2o(C, V-iio) + 2i?2o(C, M + 3i?3o(C C, 0,02. (5.6) 

The quantities in these expressions are defined as follows: is the adjoint eigenvector 
to C, the operators Rik are given by, see [H p. 95-96], 

Roi :=a^,i?(0;0) =0, 

R2o[v,w] ^D^R{0;0)[v,w] 

_ _1 /'A-^{vi+ V2) ■ Vwi + A-^{wi + W2) ■ Vwi\ 

2\A-^{vi+v2)-\7w2 + A-^lwi+W2)-Vv2)' (5.7) 



V^W2 



R30 = ^D^R = 0, 



and the functions 'i/'ijfci from the expansion of ip, are the unique solutions to 

(2c.-i-L)^2oo = i?2o(C,C), (5.8) 
-LVJiio = 2i?2o(C,C)- 

Computation of a. Since Rqi — and ker(L) = {0}, — L-iAooi = ^01 implies 
V'ooi — 0. For this result the parameter fi2 is held fixed at zero so that, using (|5.5p . 
the coefficient a of the reduced system (j5.2p is 

a= (i?n(C),C)2 = ^((51,1,0)^,02, (5.9) 
^2 

where C* is the adjoint eigenfunction, satisfying 

L*C --ic^C , (C,02 = l. (5.10) 
with the adjoint operator of L given by (using integration by parts) 

L w = V ^9 ] , V e Y2 

\ -^BV + V\I^V2 ) ' 

^''^'''^^^ ^ %L,.2 (^i-^2,gfc)gfc(a:). 
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The critical adjoint eigenfunction C*, as any eigenfunction of L*, has the form = 
T]g„i{x), where 77 = G is an eigenvector of M*^ derived from (j4.8p . If 

m 7^ (1, 1), then {(, C*)2 = 0, therefore m = (1, 1), and hence ^7] = —iujrj so that 
from (511,(711) = L1L2/2 and (|5.9p we infer 

a = < V^i*- (5-11) 
Due to dHH), there is ^ € such that 

A'h,,^^iujt^ = {e,ef, (5.12) 

and using (|5.ip at the bifurcation points AT = j = 1,2, we have 

c2 =ci(2c2-ci). (5.13) 
Together with equation (|5.10p we readily check that 

/,r • /cii-C2i-C3 -C2i \^ „ 

^'''^^^-^"^^^l C2i -Cii + C2i-C3j^ = « 

^ ^'^ ^ V C2I Cil - C2I- C3J 

2 

L1L2 

Due to (|5.13l) . the eigenvectors can be chosen as 



C11-C21-C3J \c1i-c2i-c3 
^^J^ 1 ^ 

L1L2 (cii - C3)(cii - 2c2i - C3) ' 



(5.15) 



where 5^0 provides the normalization. Therefore, (|5.1ip yields a = irc^LiSi as 
claimed. 

Computation of b. We first show ■0200 = 0; recall (|5.8p . Thanks to (|5.12p . ({x) = 
£,gi,i{x) and for fc G N* x Z we have 



L2 ki 



A,g,ix) = —l^-^^^M^). (5.16) 



A direct calculation yields i?2o(CiC) = 0. Since ker(2ia; — L) = {0} on account of 

Theorem 14.41 the equation for 7/1200 h'om (|5.8p implies V200 = 0. Together with 
i?30 = and (j5.6p . this means 

fe=(2i?2o(C,^iio),02. (5.17) 

Next, we compute Viio using (|5.8p . From = ^gu and (|3.4p . (|4.6p . (|4.5p as well as 
(|5.16p . straightforward calculations give 

- 2i - 

-L?/>iio = 2i?2o(C,C) = 17^:^^52,0. 
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Since the eigenvectors ((?fc)fe6N. xz of L are mutually orthogonal and M2,o is a multiple 
of the identity, we have that i/'iio = ct^52,o, where 



Li L2 



It follows, after straightforward calculations, that -R2o(C7 V'lio) = /^^i, 1^2,0? where 

/3 = " T.„ , AT., 



Substitution into (|5.17p yields 



b = (2i?2o(C,'0iio),C*>2 = 2/3 • 77(51, 102,0, 51, i> = --y^/S • t?- 
Finally, we use that ^ • 77 = 0, see (|5.15p . and together with £, -rj = xTlI obtain 



which concludes the proof. □ 

We now turn to the proof of CoroUarv 15.21 

Proof. [CoroUarv 15. 2] From (j5.3p and (j5.ip we readily check 6 < 0. 

Writing (j5.3p in terms of Cj and using (|5.15p . a straightforward calculation gives 

5i(a) = . ■ ^(C2 - ci). (5.18) 

^2 |(Cl - C3l)(Cil - 2C21 - C3)\ 

Thanks to Cj > 0, j = 1, 2, 3, all factors in this expression are positive, except possibly 
the last one, and therefore the sign of ^{a) is the sign of C2 — ci. Note that 

^{c2~ci) = dATd{AT,l)/Ll (5.19) 
^2 

and that the quadratic polynomial d{-, 1) has negative quadratic coefficient. Therefore, 
ci < C2 at AT = ATi and so 3fi(a) > 0, while at AT = Ar2 we have C2 < ci, hence 
3?(a) < 0. We readily compute that S(a) = cl^{a) / {03(02 - ci)) = 0{i^^). 

In conclusion, there are generic supercritical Andronov-Hopf bifurcations at both 
cndpoints of the instability region. As usual, the local invariance of the center manifold 
from Theorem 15.11 implies the existence of the claimed heteroclinic orbit between the 
unstable steady-state and the stable limit cycle, contained in the center manifold. 

Now consider the behaviour of a and b for small viscosity < ^ 1. With 
04 := ^ we get C2, C4 = 0(1), and 

2ttc?. 

a = 



b = - 



L2{oi + vc4\){{ci - 2c2)i - voi) 
Ll oj + v^ol 



2.2 (5.20) 
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Left endpoint of the instability region: ATi. Inspecting the formula for (IaT we 
find ci = C2 — \/cl — v'^c\, where c| — v'^c'l > by (|4.13|) . Hence, 



ZC2 



and wc obtain 



4„2 

4 



Therefore the radius of the stable limit cycle \z{t)\ for sufficiently small is 

2. ^(4^?+^i)c.Vy^^(,,). 



Right endpoint of the instability region: AT2. Here ci = C2 + ycf — i^^^, therefore 
ci = 2c2 + ©(i^^) and so 



c|L2 f TT^{ALl + Ll) V 

hence the radius of the stable limit cycle for small — /ii is 



4Lf + Li 



This concludes the proof. □ 



We finally provide the proof of Theorem [ 

Proof. [Theorem I5.3j In order to unfold in we cannot cite a center manifold 
theorem from [3] verbatim. The reason is that /i2 modifies the second order derivative 
terms, but the results in [4] are formulated only for parameter dependence of lower 
order terms. However, as pointed out in [4l Remark 3.7], there is no problem, if 
the domain of L is independent of the parameter. This is the case here as long as 
V = i'crit(l) ~ M2 > 0; which is valid for the purpose of unfolding from v = Vcnt- 
More precisely, the proof of [H Theorem 3.3, p. 46], which considers the phase space 
extended by the unfolding parameter space, applies as follows for i'crit(l) > /^i due 
to the linearity in /ij. Set /i = (/ii,/i2), u = (m, /i) and Lu = (L + /.tiO^^^ 0)"^ ~ 
^\V^u, 0) as well as R{u) ~ {R{u), 0). (We use ^2 the parameter instead of ^2 only 
for more pleasant reduced equations.) For the extended problem, the parameter-free 
center manifold theorem applies [4j Theorem 2.9]. 

Therefore, as in the first part of the proof of Theorem 14. 4[ we obtain existence 
of the center manifold and the coefficient b is unchanged. Let A denote the real 
coefficient of z in the vector field on the center manifold. It remains to derive the 
claimed a j -dependent form 

A aini (02/^2 - flsA^i) + 0{nl + \^J.l^J-l\). 

For this we simply note that in the present case, (|5.9p is replaced by the more general 
form 



A - (i?n(C)Ai,02 = ma - M2(V'C,C)2, 
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where C = ^51,1. Using V^gi^i = — tt^ + -^j gi^i as well as (C, C*)2 = 1, we obtain 



1 4 



A = ^iia + iil^'' ( T2 + 72 ) = Mia + ©(/i^), (5.21) 



with a from Theorem 15.11 whose dependence on /12 is considered next. Recall that 



/^2, with /Z2 = giving equality in (|4.13p . Hence, 



2 



I'crit 



7rLiL2C§ ' 
where £3 stems from writing 

d(Ar, l)/L\ = (2Si - £2Ar)Ar - 1^252^ 

with suitably defined Cj, 7 = 1, 2, 3 (note the relation to Cj in (IS.ip V Then d{AT, 1) 
gives 



ATcrit = Ci + yAi2\/2l'crit " m1 

Using dnUl), ((5l9l) with AT = ATcHt + Mi then yields 



a = fli ^£3/12-^/21^^-^1 - 2c2Mi 

The above formula for Vcnt and expansion in fi2 ~ gives claimed form of A, when 
substituting the resulting a into (j5.2ip . 

The bifurcation scenario can be immediately read off the reduced vector field. □ 

6. Travelling wave bifurcation. As mentioned in the introduction, due to the 
translation symmetry in X2 , the Andronov-Hopf bifurcations correspond to periodic 
travelling wave bifurcations. Specifically, each periodic orbit is a steady state in a 
comoving frame 1/2 ~ X2 — st for certain s. While this is somewhat folklore, for 
completeness we give some details. The converse is clear: periodic travelling wave 
bifurcations imply Andronov-Hopf bifurcations. 

First note that the effect of the co-moving variable is the introduction of an ad- 
vection term sdy^ on the right hand side of the first two equations in ()l.ip . Therefore, 
the linearization Mk is replaced by 

Mk,s ^ Mk + sCi{k)ld, 

where Ci(fc) = 2iT\k2/L2- Hence, if \k is an eigenvalue of then Afe -I- sCi(fc) is an 
eigenvalue of Mk_s and choosing critical fc2 = +1, the frequency at bifurcation w is 
replaced by w + s27r/L2- The reduced equation on the center manifold then reads 

z = i(a; + s2ti/ L2) + /iia + 6z|zp, 

where a and b are unmodified since the matrices made of Cj in (|5.14p do not depend 
on s. Hence, for s = s.^ ;= —ljL2/2'k we find steady state supercritical pitchfork 
bifurcations. Note the choice fc2 = —1 reverses the sign of simply leading to the 
complex conjugate equation. 

This argument is slightly incomplete since the spectrum of the modified L pos- 
sesses a double zero eigenvalue at s = . Hence, the coefficients on the center manifold 
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are not immediately given by the Andronov-Hopf case used above. However, the re- 
duced vector field on the 2D center manifold of the double zero eigenvalue reduces to 
a scalar equation, undergoing a pitchfork bifurcation, precisely due to the translation 
symmetry. In polar coordinates of the Hopf normal form, this is due to detuning 
the trivial angular equation, co- rotation with velocity sCi(l), to stationarity. Such 
reductions due to continuous symmetry also hold in more abstract contexts, see, e.g., 
Theorem 2.18 of [4], where an additional rcfiection symmetry is assumed. 

In the context of travelling waves, let us briefly take the perspective of pattern 
formation, for which the infinite strip x € [0, Li] x R is the natural domain here. The 
linear stability analysis of the laminar in this case involves the eigenvalues from 
21 with continuous and rescaled these are eigcnmodes in the essential spectrum 
given by with k ~ (fci, ^2^:2), (fci, fc2) G x M. In particular, the critical modes 
can only be A^^ L2k2)^ ^2 G I^- 

We are then lead to search for pattern-forming instabilities, and indeed, the sys- 
tem easily allows for the analogue of Turing- Hopf instabilities from reaction-diffusion 
systems, which is also well known in fiuid dynamics, for instance Rayleigh-Benard 
convection. A detailed analysis is tedious, and we only give a numerical example in 
Figure 16.11 which is derived from that in Figure I4.2f a') . Here the critical modes at 
onset of the instability on the infinite strip have wavenumber near 0.75. The periodic 
solutions of SjSj alias, wavetrains, are a signature of the bifurcating continuum of peri- 
odic solutions. The fact that these are supercritical suggests supercritical Turing-Hopf 
bifurcations. 

Re A 

0.02 h 



-0.02 

-0.04 
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Fig. 6.1. A Turing-Hopf instability for parameters as in Figure \4-2}[ a): real parts of eigenvalues 
A^, k = (1, fc2), as functions of k2 for AT = 0.05 (stable), AT = 0.08 (near bifurcation) and 
AT = 0.2 (unstable). 

7. Nonlinear Instability. In this short paragraph, we give some details on 
the fact that the linear instability of the laminar state pss is indeed an instability 
for the nonlinear equation uniformly in ly. Roughly speaking, this means that in 
(ATi, AT2), there are initial data which are arbitrarily close to the steady state and 
which get "far" from it exponentially quickly. We thus assume ATi < AT2 and take 
AT G (ATi, AT2). This means that cr+(L) ^ (sec the proof of Theorem |13D. 

In the parabolic formulation p.Sp , the sectoriality of L allows to apply the well- 
known nonlinear instability results from [5] for spectrum in the right half plane. How- 
ever, this heavily relies on v > and the following docs not. Furthermore, the result 
given for the specific case here is actually stronger than the general ones in [S]. 

As in [21 Theorem 6.1], the following instability result holds for > 0: 
Theorem 7.1. Suppose AT € (ATi, Ar2)- There exist constants (5o,?yi,?y2 > 
such that for any < S < 5q and any s > there exists a solution (p* , E) to (|l.ip 
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with 11/3(0) — Pss\\h^ < <^ but such that: 

Wpits) - PssIIl^ > Vi and j|£'(t5)||L2 > t]2, 

with ts = 0{\logS\). 

Since the proof is almost identical to that of Theorem 6.1 in [3], we refer to this 
paper for details. The main idea is to apply the method of Grenier [5]. 

8. Global Nonlinear Stability. Let us now investigate the stability of the 
steady state pss, outside of [ATi, Ar2]. Wc first state the results and then give the 
proofs. 

The key point is the following energy identity: 



Lemma 8.1. For any initial data po G L° 
the solution of the system (|l.ip 



we have the following estimate for 



£{t) \\p-pss\\h- 



LiAT 
t 

2v 



V^VVdx 







-2 



LiAT 
< £{0) = ilPo 



P 



ds 



pssWh 



LiAT 



SS/IIL2 



l^V^t^ol^dx, 

lV(p+-p+)ll 



L2 



wtth Hp - pssWl^ = \\p+- ptsWh + \\P~- P7s\\h, II v(p - p 

l|V(p" - p-)\\h c^nd VV = V(V2)"1(P+ + - l). 

Remark 7. We take the opportunity to point out an error in the energy of 

Theorem 5.1 J: in equations (5.1) and (5.2) of this paper, there is a factor 2 which is 

missing in front of \\'V\'^dx. 

We shall use in the following Poincare type inequalities: 

Lemma 8.2. With the same notations as before, we have, for any t > 0: 



9 r 2 

\yv\\l. < '-^ 



\\P- PssWl^, 



< ^\\Vip-Pss)\\t2. 



(8.2) 
(8.3) 



As a consequence of the energy identity, we can prove L^-return to equilibrium, with 
exponential (and explicit) speed, for negative or large enough AT. 

Theorem 8.3. If AT < or AT > AT* := then the steady-state pss is 

globally asymptotically stable in L^, with exponential convergence given by — 27r^y^. 

Remark 8. Notably, the constants in all these results are independent of L2, 
and the convergence rate is larger on thinner domains (with smaller Li), but also 
balancing with viscosity. 

Recall that by Lemma \4.3[ for fixed £ ~ L2/ Li the higher instability threshold AT2 
satisfies AT2 < AT^ = but that limf_>._|-oo AT2 ~ AT* if v ~ o(£^^). Hence, the 
global threshold AT* is also linked to linear instability. 



Let us now Lemma [8. II Lemma 18^ and Theorem 18.31 

Proof. [Lemma l8.1j The proof follows from computations that are similar to those 
that can be found in [3], for the model without viscosity (that is = 0). We keep the 
notations of Section [31 
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Taking the scalar product with u := (ui,it2) in the transport equations satisfied 
by ui and U2 in (|3.ip . and integrating with respect to x entails: 



Id 

7i 

2dt" ' 



L2 = /" — -ui dx — f — -U2dx + i' ( ui^l'^uidx + v ( U2'^'^U2dx. (8.4) 
Jo, Jn Li Jq 



Note indeed that due the periodicity with respect to X2, the following contribution 
vanishes: 

/ dx2UiUidx — j -dx^uldx = = j dx2U2U2dx. 
Jn Jn 2 Jn 

Likewise, with Green's Formula, using div E-^ = and E2 = —5^2^ = on xi = 0, Li, 
we have (for i = 1,2): 

/ E-^ ■ Wu, u,dx = - / E-^ ■ V{uifdx = 0. 
Recall an identity proved in [3l Lemma 5.1]: for any t > 0, there holds 



E2Uidx = E2U2dx. (8-5) 
n Jn 

For the sake of completeness, we quickly reproduce the proof. Observe that 

/ E2{u2~u\)dx^ I E2 {ui + U2 — 2ui) dx 
Jn Jn 

= [ E2 {~V^V - 2ui) dx (8.6) 
Jn 

= —2 / E2Uidx. 
Jn 

Indeed, relying on the periodicity in the X2 direction and since S^^jT^ = on xi = 0, i, 
we get: 

/ dx^V^'^Vdx = - [ dx^'^V ■ \7Vdx + [ div{dx2V\7V)dx 
Jn Jn Jn 



=0 



( ]dx = 0. 



This completes the proof of (|8.5p . Therefore, we have: 



^2 f ^2 r, f ^2 

ui dx + — U2 dx = —2 / — ui dx. 

Li Jo Li Li 
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Now compute, using the equations satisfied by (ui, U2,V): 

/ E2U1 dx ~ / Vdx2Uidx — / div{Vuie2)dx 
Jn Jn Jn 



(at(ui + U2) + • V(ui + U2)) dx 

+ ^ ^ -r-y9.,U2da; " ^ VV^im + U2)dx 

= ^ / (StV^y - ■ V(VV)) dx-^ j VV\u^ + ^i2)da;. 

(8.7) 

Observe that by Green's formula: 

/ -V {dt^^V - ■ ^{^^V)) dx^~( f \^V\''dx\ . 
Finally, using again (|8.5p . we have 



E2 J f E2 , 
■ui dx — —U2 dx 



1 d f f ,„,,,, \ . 2v 



Li{T+ -T-)dt 



Note that using Green's formula and the Poisson equation satisfied by V , we have the 
identities: 

/ V\J'^ {Ui + U2)dx = / \J'^V{UI+U2)dx ^ -\\ui+U2\\l2. 

Jn Jn 

and 

V i uiV'^uidx + v / U2V'^U2 dx = —ly / \\7ui\'^ dx — v / |Vw2|^dx 
Jn Jn Jn Jn 

Gathering all pieces together, we have proved that ■^S{t) = 0. □ 
Let us now prove the Poincare inequalities of Lemma 18.21 

Proof. [Lemma [HH] We only prove (|8.2D f (|8.3|) can be treated similarly). Using 
the orthogonal basis (|4.3p . we write: 

Ui+U2= ^ aki^kidk- 

fci6N,,fc2eZ 

Recall the Poisson equation satisfied by V: 



'VV = Ui + U2. 



This yields: 



2 / i£i 
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Therefore, 



/ iVFpdx < -^\\u, + u,\\l < -^i\\u,\\l, + \\U2\\1), (8.8) 

which proves (|8.2p . □ 

Gathering all pieces together, we can now prove Theorem 18. 31 

Proof. [Theorem 18 .3] Using the energy identity (|8.ip and applying the Poincare 

inequality (|8.2p we get: 



Wp-pssWh < \\piO) - PssWh + 



LiAT 



- 2v 



-4 



LiAT 



\\P- Pss\\h + \\V{p- Pss)\\l2 



ds 



1 4^2 



- 2v 



Hence, using the Poincare inequality 



t r A 



LAT 



\p-PsA\t^ + \\'^{p-Pss)\\t2 



ds. 



4il \ 11 .,2 



( ' - 


-) 




\LiAT 




L 



< \\p{0)~Pss\\l2+2l 



As a consequence, by Gronwall inequality, wc obtain L^-stability and L^-rcturn to 
equilibrium, provided that 



AT < 0, or AT > 



4Ti 



which in particular implies — ^^^^^ > 0. More specifically, we have: 

\\p - PssWh < \\p{0) - PssWh cxp (-7t) , (8.9) 



with 7 := 2u - ^) (f 



4Li \ 
7t2AT} 



2ufy > 0. □ 
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